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I.  INTRODUCTION  AND  SUMMARY 


The  Department  of  Defense  (DoD)  is  well  aware  of  the 
high  costs  of  "scientific"  software.  The  trend  of  increasing 
software  costs  is  well  established,  and  for  most  computer  systems 
the  software  accounts  for  50-90%  of  the  total  system  cost.  In 
the  case  of  hand-held  calculators,  the  situation  is  even  more 
extreme  as  the  hardware  costs  are  usually  well  under  $500  cind 
software  costs  may  reach  into  the  tens  of  thousands  of  dollars 
per  program.  As  with  most  computers,  these  programs  are  machine 
dependent,  that  is,  they  will  only  run  on  the  machine  for 
which  they  were  designed.  The  rapid  developments  in  the  semi¬ 
conductor  industry  has  created  a  situation  where  new  calculators 
are  being  introduced  every  1-3  years.  It  is  desirable  to  keep 
these  programs  running  on  the  current  state-of-the-art  machines, 
both  from  the  standpoint  of  decreased  execution  time,  but  also 
parts,  service  and  technical  support  from  the  manufacturer.  Be¬ 
cause  these  programs  are  machine  dependent,  a  substantial  amount 
of  code  conversion  is  required  to  keep  these  programs  operational 
on  the  most  current  machines.  This  conversion,  if  done  by  hand, 
is  costly,  error  prone  and  time  consuming. 

The  Minicomputer  Assisted  Reprogramming  System  (MARS)  was 
conceived  in  1977  to  address  this  problem.  MARS  provides  a  hedge 
against  rising  software  costs  and  early  program  obsolescence  by 
providing  machine  assisted  software  translation  of  calculator 
programs  to  the  new  calculator  languages  rather  than  a  complete 
recoding  of  the  program.  In  addition  to  helping  lower  the  costs 
of  implementing  existing  programs  on  new  machines,  MARS  can  pro¬ 
vide  a  quicker  turn-around  time  than  recoding  and  provides  a 
means  to  pool  software  resources  for  currently  available  calcu¬ 
lators.  Prior  to  the  second  year  of  funding,  the  MARS  concept 
was  reviewed  in  light  of  total  DNA  software  requirements.  The 
results  of  this  review  revealed  that  the  scope  of  the  project 
should  be  broadened  to  include  both  general  and  specific  software 


support  tools.  As  with  standard  computer  software  development 
efforts,  software  tools  are  one  way  to  help  control  costs,  de¬ 
crease  development  time  and  insure  a  quality  product. 

Prior  to  the  MARS  project,  no  development  or  support  tools 
for  calculator  software  were  available.  We  undertook  the  design 
and  development  of  these  tools  to  assist  the  calculator  program¬ 
mers  in  their  day-to-day  tasks.  These  tools  include  an  automated 
"time- space"  tradeoff  code  optimizer,  TI-59  and  Constant  Read 
Only  Memory  (CROM)  simulator,  an  algorithm  development  facility 
that  consists  of  a  sophisticated  curve  fitting  package  and  sup¬ 
port  routines  and  a  graphics  package  that  provides  for  curve  fit 
verification  and  assists  in  program  verification. 
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II.  HARDWARE  AND  SOFTWARE 


Early  during  the  first  year  of  funding  of  MARS  all  of  the 
available  options  for  providing  computational  services  were  re¬ 
viewed.  This  included  a  survey  of  timesharing  services  and  a 
study  of  available  low-cost  minicomputer  based  systems.  Acqui¬ 
sition  of  a  minicomputer  was  proposed.  The  reasons  included 
higher  system  availability,  hands-on  management  of  off-line 
storage,  the  ability  to  interface  additional  peripherals,  and 
if  the  requirements  for  the  handling  of  classified  data  arose, 
physical  accountability  of  the  computer  data.  Appropriate  candi¬ 
date  computer  systems  were  reviewed,  bids  solicited  and  a  system 
was  chosen.  The  MARS  translation  system  is  based  on  the  Digital 
Equipment  Corp.  PDP11T03.  The  11T03  system,  as  ordered,  consists 
of  the  LSI-11  CPU;  KEV-11  bootstrap  and  diagnostic  ROM;  56  kb  of 
MOS  RAM;  2  RK05  disk  drives  and  controller  organized  as  one  re¬ 
movable  disk  and  one  fixed  disk  in  three  logical  platters  for  a 
total  of  7.5  Mb  of  storage  on-line;  2  DLV-11  asynchronous  serial 
interfaces;  one  VT-52  CRT/keyboard  and  one  LS-120  printer/keyboard. 

During  the  second  year  of  funding  additional  hardware  was 
added  to  the  11T03  to  provide  graphics  input  and  output  cap¬ 
ability.  This  hardware  consists  of  a  Summagraphics  BP-11  "Bitpad" 
digitizing  tablet,  cursor,  power  supply  and  RS232  interface;  a 
DLV-llF  single  channel  RS232  interface;  Hewlett  Packard  7225A 
X-Y  plotter  and  IBV-llA  IEEE-488  interface.  The  11T03  system 
with  peripherals  has  functioned  normally  and  has  proven  to  be 
quite  reliable. 

All  MARS  software  is  written  and  designed  for  interaction 
with  DEC'S  RT-11  operating  system.  Version  3-02  is  used  by  MARS 
and  provides  all  functions  normally  associated  with  operating 
systems  including  overlay  capability  and  disk  file  management. 

From  the  outset  of  the  MARS  project  the  requirements  for  port¬ 
ability  and  maintainability  of  the  applications  software  have 
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been  undefined.  As  a  consequence,  all  routines,  with  a  few 
exceptions,  have  been  coded  using  DEC'S  version  of  ANSI  standard 
Fortran  IV  for  portability,  readability  and  maintainability. 

Initial  MARS  software  development  was  done  using  DEC's 
Fortran  IV  V2.04,  but  the  discovery  of  several  compiler  "bugs" 
and  shortcomings  has  forced  us  to  upgrade  to  version  2.1.  In 
addition  to  Fortran,  the  MARS  algorithm  verification  facility 
uses  DEC'S  single  user  BASIC  version  1.  All  MARS  files  are 
written  in  the  Standard  RT-11  file  format  and,  where  appropriate, 
are  written  using  ASCII  code  for  human  readability. 
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III.  TECHNICAL  APPROACH 


MARS  was  originally  conceived  as  a  language  translation  aid 
to  help  provide  a  hedge  against  the  obsolescence  of  calculator 
based  scientific  software.  As  the  system  design  progressed 
and  more  experience  was  gained  in  the  processing  of  calculator 
software,  it  became  obvious  that  the  original  concept  could,  and 
should,  be  expanded  to  provide  automated  support  in  other  areas 
of  the  calculator  software  development  task.  During  the  second 
year  of  funding,  additional  features  were  incorporated  into  the 
existing  software  system.  As  can  be  seen  by  Figure  1,  the  MARS 
system  currently  consists  of  three  functional  areas:  the  Lan¬ 
guage  Processing  aids;  Data  Processing  and  Algorithm  development 
aids;  and  Program  Development,  Verification  and  Support  aids. 

A  description  of  the  purpose,  implementation  and  a  brief  functional 
description  of  each  of  the  component  modules  follows.  Appendices 
A-1  through  C-4  contain  listings  of  these  modules  and  a  more  com¬ 
plete  description  of  the  operation  of  the  module  as  program  com¬ 
ments.  Where  appropriate,  examples  will  be  included  in  this 
section  of  the  report  to  demonstrate  critical  technical  points. 

III-A.  MARS  Language  Processing  Modules 

This  group  of  Fortran  IV  modules  comprises  the  bulk  of  the 
implementation  of  the  original  MARS  concept.  The  basic  idea  was 
to  provide  an  expedient,  machine  assisted  method  for  making 
existing  calculator  programs  executable  on  newer  or  more  sophis¬ 
ticated  calculators.  At  the  outset  of  the  project,  only  two 
manufacturers,  Hewlett  Packard  (HP)  and  Texas  Instruments  (TI) , 
were  producing  hand-held  personal  programmable  calculators  (HHPPC) 
with  enough  capabilities  for  DNA's  software  tasks.  Our  studies 
revealed  vast,  fundamental  differences  in  the  programming  and 
operation  of  these  two  classes  of  machines.  Furthermore,  it 
was  not  clear  whether  other  manufacturers  might  market  models 
with  the  necessary  capabilities.  As  a  consequence,  it  became 
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MARS  CAPABILITIES  ORGANIZED  BY  FUNCTIONAL  AREAS 


Figure  1. 
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apparent  that  our  solution  to  the  language  translation  problem 
would  have  to  be  general  purpose  in  nature  and  very  expandable. 
Additional  system  design  parameters  were  gained  through  the  in¬ 
sight  that  if  one  were  to  implement  a  single  purpose  "point-to- 
point"  translator  to  provide  the  ability  to  go  from  machine  A  to 

machine  B,  and  there  exist  N  machines,  then  the  number  of  trans- 

2 

lator  modules  required  increases  as  N  -N  (Figure  2) .  However, 
if  one  could  design  a  general  purpose  intermediate  language,  all 
translations  could  be  made  through  this  language  and  the  number 
of  translator  modules  only  increases  as  2N  (Figure  3) .  During 
the  analysis  phase  of  MARS,  HP  and  TI  were  each  releasing  new 
machines  on  a  two-year  design  cycle,  one  year  out  of  phase  with 
each  other.  The  result  was  that  a  new,  more  powerful  machine 
was  being  released  each  year  from  alternating  manufacturers. 

This  situation,  coupled  with  the  fact  that  each  manufacturer  had 
two  machines  in  the  field,  pointed  strongly  to  the  intermediate 
language  solution  being  the  best.  As  an  additional  benefit  to 
this  solution,  the  intermediate  language  provides  a  convenient 
storage  mechanism  for  translated  programs  and  provides  the  ability 
to  modify  and  optimize  the  program  in  its  intermediate  form. 
Furthermore,  this  solution  provides  for  unlimited  future  expansion 
in  terms  of  the  number  of  source  and  target  machines  that  can  be 
handled. 

Several  existing  languages  were  considered  as  candidates 
for  the  intermediate  language.  All  were  discarded  due  to  data 
formatting  problems  and  overall  complexity  of  the  finished  program. 
We,  instead,  opted  to  design  our  own  intermediate  language.  Our 
prime  design  considerations  were  human  readability,  simplicity 
and  ease  of  translation.  The  completed  design  of  SIMPLE  (for 
Special  Interpretative  Mathematical  Programming  Language  Experi¬ 
ment)  contains  constructs  for  all  calculator  based  mathematical 
operators,  addressing  modes,  and  transfers  of  address.  The  most 
complex  expression  allowed  in  SIMPLE  is  of  the  form  "A  =  B  +  C". 
Expressions  that  are  more  complex  are  decomposed  by  parsing  until 
they  are  in  this  elementary  form.  The  main  advantage  of  this 
form  is  that  it  promotes  human  readability  and  facilitates  machine 
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Point-to-Point  Translation  Approach 
Figure  2. 
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Intermediate  Language  Translation  Approach 


Figure  3. 


manipulations  of  expressions.  Generally  speaking,  the  trans¬ 
lation  process  requires  two  steps  to  complete  a  source  to  target 
translation.  Step  one  is  to  translate  from  the  source  code  syn¬ 
tax  into  a  form  compatible  with  the  syntax  of  SIMPLE.  The  second 
step  is  to  translate  from  SIMPLE  to  a  form  acceptable  to  the  tar¬ 
get  machine.  This  arrangement  provides  for  stand-alone  operation 
of  each  translation  module.  All  processing  modules  communicate 
with  each  other  via  disk  files  using  the  standard  MARS  file  for¬ 
mat.  Since  each  module  is  designed  for  stand-alone  operation, 
and  no  two  modules  are  core  resident  simultaneously,  an  unlimited 
number  of  modules  can  be  added  between  any  two  processing  end 
points.  This  concept  insures  expandability  to  any  desired  degree, 
and  will  help  protect  the  system  from  obsolescence. 

The  language  processing  portion  of  MARS  has  been  divided 
into  four  modules.  The  modules  provide  the  ability  to  move  source 
code  to  and  from  the  SIMPLE  language.  The  modules  provide  the 
capability  to  translate  from  HP67/97  code  to  SIMPLE  (HP67/97 
decompiler);  from  TI59  code  to  an  intermediate  in-line  TI59  code 
(TI  preprocessor) ;  from  TI59  in-line  code  to  SIMPLE  (TI59  de¬ 
compiler)  ;  and  from  SIMPLE  to  TI59  (TI59  recompiler) .  Due  to 
the  radically  different  methods  of  problem  representation  of  the 
two  manufacturers,  i.e. ,  RPN  vs.  AOS,  large  difficulties  in  the 
translation  process  have  been  encountered.  Both  systems  perform 
problem  evaluation  in  an  interpretive  manner  using  stack  oriented 
processing  for  expression  parsing  and  the  actual  numerical  evalu¬ 
ation.  In  the  HP  case,  the  stack  processing  in  an  integral  part 
of  the  machine  and  actively  manipulated  by  the  programmer.  All 
operators  are  executed  as  they  are  encountered  and  operate  on 
data  that  resides  in  the  stack  at  the  time  the  operator  is  encoun¬ 
tered  in  the  program.  For  the  TI  machines,  stacks  also  exist, 
but  are  implicit  in  the  operation  of  the  machine.  Because  the  AOS 
representation  of  the  task  is  not  immediately  executable  by  the 
machine,  a  large  amount  of  parsing  must  be  done  to  allow  evalua¬ 
tion  of  the  expression.  This  parsing  revolves  around  the  use  of 
two  stacks,  one  for  data  and  one  for  pending  operations  (in  AOS 
not  all  operators  are  executed  as  they  are  encountered. )  The 
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problem  centers  around  the  fact  that  both  machines  allow  trans¬ 
fers  of  address  at  any  place  in  the  code  —  even  within  the  middle 
of  equations  or  expressions.  Although  this  is  a  very  powerful 
tool  when  used  properly,  it  presents  a  formidable  task  for  the 
language  translators.  The  evaluation  of  any  expression  is  deter¬ 
mined  by  the  "data"  in  the  stack (s),  and  since  the  translators 
must  perform  a  partial  evaluation,  the  translation  process  is 
determined,  or  at  least  affected,  by  the  data  in  the  stack (s). 

The  problem  is  further  compounded  by  the  fact  that  TI  allows 
transfers  of  address  with  pending  data  and  operators.  The 
straight-forward  solution  is  obvious  but  results  in  an  unaccept¬ 
ably  large  translated  program.  This  method  calls  for  the  evalu¬ 
ation  of  each  alternate  course  of  program  flow  in  the  code, 
taking  into  account  the  possible  changes  in  the  values  in  the 
stack (s).  Each  alternate  course  (branch,  jump,  subroutine  call 
or  reset)  requires  evaluation  all  the  way  to  the  completion  of 
the  program.  The  result  is  that  a  relatively  straight-forward 
input  code  is  transformed  into  a  huge  "tree"  with  a  "branch" 
"growing"  to  the  end  of  the  program  at  each  alternate  course  of 
action.  For  the  HP67/97  the  problem  has  been  solved  in  an  expe¬ 
ditious  manner.  It  involves  five  complete  "translations"  and  a 
substantial  amount  of  code  compaction,  including  both  local  and 
global  variable  minimizations.  With  the  exception  of  testing  all 
possible  code  sequences,  the  HP67/97  decompiler  is  complete. 
Furthermore,  our  tests  show  that  for  many  code  sequences,  the 
translated  code  is  optimal  (i.e. ,  no  unnecessary  steps) .  The 
TI59  language  processors  are  near  complete.  With  the  exception 
of  extensive  testing,  the  TI59  decompiler  is  complete.  However, 
this  presupposes  the  existence  of  in-line  code.  This  conversion 
to  in-line  code  is  the  responsibility  of  the  preprocessor.  Most 
of  the  services  required  of  the  preprocessor,  i.e.,  input  editing 
constant  and  digit  manipulation  and  subroutine  expansion,  are 
available.  However,  due  to  the  extreme  complexity  of  the  "pending 
op"  problem,  this  feature  is  unavailable.  Furthermore,  since  the 
need  to  translate  large  numbers  of  TI59  codes  to  other  machines 
has  not  yet  occured,  this  task  has  been  deferred  and  resources 
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have  been  redirected  to  other  more  pressing  problems.  See 
Appendix  A2  and  A3  for  more  details  on  the  TI59  preprocessor  and 
decompiler. 

The  last  segment  of  the  language  processing  modules  is  the 
TI59  recompiler.  The  recompiler  takes  as  input  the  intermediate 
language  representation  of  the  program  and  converts  it  to  a  form 
acceptable  to  the  TI59.  The  recompiler  attempts  to  exploit  all 
the  capabilities  of  the  TI59  and  performs  an  "inverse  parsing" 
function  in  that  it  recreates  statements  with  parentheses  under 
the  precedence  constraints  of  the  TI59.  The  output  of  the  re¬ 
compiler  is  an  unoptimized,  linear  program  sequence  that  is 
executable  but  inefficient.  This  module  is  complete,  except  for 
two  argument  numerical  functions  and  indirect  transfers  of 
address.  It  is  felt  that  these  problems  are  best  resolved  by 
human  assistance,  since  complete  automation  of  these  functions 
would  be  excessively  expensive.  The  recompiler  has  been  tested 
for  elementary  code  sequences  required  for  validation  of  the 
algorithm.  Due  to  the  excessive  amount  of  time  required  to 
synthesize  test  programs,  evaluate  the  source  and  output  for 
validity  and  efficiency,  it  is  recommended  that  further  work  on 
this  phase  be  deferred  until  such  time  as  "production"  trans¬ 
lation  programs  become  availcible.  This  method  will  save  the 
cost  of  test  program  synthesis  and  provide  a  more  realistic  basis 
for  validation. 

III-B.  MARS  Data  Processing  and  Algorithm  Development  Modules 

During  the  second  year  of  funding,  the  scope  of  the  MARS 
project  was  expanded  to  include  NWE  programming  aids.  More 
specifically,  the  four  modules  that  comprise  this  functional  group 
combine  to  provide  the  NWE  programmer  with  a  comprehensive  set  of 
tools  to  assist  him  in  numerical  modeling  of  nuclear  phenomenology. 
Many  of  the  tasks  encountered  by  the  NWE  staff  entail  the  repre¬ 
sentation  of  empirically  and  theoretically  produced  data  in  a 
form  that  will  fit  in  the  HHPPC.  More  often  than  not,  this 
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results  in  the  approximation  of  sets  of  graphs  and  tables  of 
numbers  by  a  family  of  equations.  This  task  is  collectively 
referred  to  as  "curve  fitting".  Manual  fitting  of  large  amounts 
of  data  is  costly  and  error  prone.  Therefore,  we  have  under¬ 
taken,  and  completed,  the  development  of  a  group  of  progreuns  to 
assist  the  NWE  staff  in  their  curve  fitting  efforts.  The  pro¬ 
cess  of  fitting  curves  can  be  functionally  divided  into  four 
parts.  The  first  task  is  acquisition  of  data  from  the  curve  to 
be  fit.  To  do  this  acquisition  manually  requires  the  use  of  a 
backlighted  table,  high  resolution  graph  paper  and  lots  of  pa¬ 
tience.  We  have  resolved  this  problem  with  the  use  of  an  elec¬ 
tronic  digitizing  tablet  and  necessary  computer  interface  cir¬ 
cuitry  and  software.  The  tablet  provides  an  11"  x  11"  active 
area  with  2200  points  resolvable  in  each  direction  resulting  in 
a  resolution  of  .0005".  Coordinate  information,  as  specified  by 
the  operator  using  a  stylus  or  cross-haired  cursor,  is  trans¬ 
formed  to  a  four  digit  ASCII  absolute  tablet  coordinate  (ATC) 
and  transmitted  to  the  11T03  via  an  RS232  interface.  Since  with¬ 
out  scaling  information  the  ATC  data  is  useless,  we  have  com¬ 
bined  several  management  functions  into  one  module  for  data  ac¬ 
quisition.  This  module,  DAM4,  for  Data  Acquisition  Module  ver¬ 
sion  four,  is  an  interactive  Fortran  IV  program  that  provides  a 
variety  of  services.  It  allows  for  calibration  of  the  tablet, 
compensation  for  translation,  rotation  and  scaling  ci  the  graph 
using  a  sophisticated  multi-point,  least  squares  coordinate 
transformation,  and  provides  for  the  entering  of  axis  label 
information  and  data  file  management  including  multiple  curves 
per  file.  The  result  of  a  session  using  DAM4  is  a  formatted 
ASCII  disk  file  containing  the  coordinate  data  for  points 
selected  by  the  user.  This  data  file  is  passed  on  to  the  next 
link  in  the  processing  chain  —  DERM. 

DERM  stands  for  Data  Edit  and  Review  Module,  and  is  an 
interactive  progreun  that  does  as  the  name  implies.  DERM  uses 
the  output  file  from  DAM4  as  input  and  provides  the  user  with 
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the  ability  to  quickly  review  the  digitized  data  for  smoothness 
and  accuracy,  perform  single  and  double  argument  immediate  trans¬ 
formations  on  the  data,  insert,  delete  and  change  data  points, 
list  data  on  the  printer  and  produce  "quick"  plots  of  the  data 
for  visual  evaluation.  Furthermore,  DERM  has  the  capability  to 
provide  for  the  linking  and  unlinking  of  related  data  files  and 
allows  for  data  inversion  on  families  of  graphical  data  in  either 
single  or  multi-file  format.  As  its  output  DERM  produces  a  data 
file  with  corrected  points  that  is  passed  on  to  the  Fitter  module. 

The  fitter  module  is  the  crux  of  the  data  processing  pack¬ 
age  in  that  this  module  is  responsible  for  the  actual  curve 
fitting.  Unlike  earlier  attempts  at  an  automated  curve  fitter, 
we  have  kept  the  human  "in  the  loop"  to  insure  a  reasonable  choice 
of  base  function  and  resultant  coefficients.  The  actual  task  of 
fitting  numerical  data  is  distributed  over  five  submodules.  The 
first  submodule  is  responsible  for  disk  file  input/output  and 
provides  the  user  with  dynamic  editing  capabilities  as  well  as 
immediate  data  transforms  and  scratch  pad  management  functions. 
Siibmodule  two  is  the  actual  curvefitting  routines.  Submodule 
three  is  the  function  form  definition  module.  The  user  has  the 
ability  to  choose  any  of  60  "prepackaged"  base  function  forms 
as  well  as  any  one  of  five  user  defined  forms.  The  functional 
form  choices  include  polynomial,  logarithmic,  normal  distribution, 
inverse  cumulative  normal  distribution  and  Legendre  polynomials. 
Submodule  four  is  the  function  evaluation  module.  In  addition 
to  being  able  to  evaluate  the  fitting  function  given  the  function 
form,  this  submodule  gives  the  user  the  ability  to  choose  coef¬ 
ficient  initial  values,  search  intervals  and  starting  values. 
Additionally,  the  search  routine  can  be  constrained  according  to 
user  specified  criterion,  like  mean  square  error  or  maximum 
error.  It  also  allows  evaluation  of  the  function  at  user-defined 
initial  values  and  intervals,  and  can  store  these  values  for  use 
by  the  plot  module.  The  fifth  submodule  is  a  data  smoother. 
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This  allows  the  user  to  represent  the  given  data  using  cubic 
spline  techniques  to  provide  for  function  evaluation  at  periodic 
intervals  of  the  independent  variables.  This,  then,  provides  an 
ordered,  periodic  sampled  data  base  for  the  fitter  submodule. 

The  final  module  in  the  Data  Processing  package  is  the  Plot 
module.  This  module  allows  the  user  to  produce  high  quality 
plots,  including  alphanvimeric  and  axis  data,  on  the  HP7225A  X-Y 
plotter.  The  plot  module  uses  as  input  files  produced  by  DAM4, 
DERM,  Fitter  or  the  user  and  generates  plots  of  the  ordered  data 
on  8*5"  X  11"  paper.  The  plot  module  attempts  to  use  the  para¬ 
meters  contained  with  the  data  file  to  define  all  the  "particulars" 
of  the  plot  —  like  size,  orientation,  axis  scaling,  etc.  The 
user  has  the  ability  to  override  the  defaults  at  any  time.  This 
module  allows  linear,  logarithmic  or  user-defined  axes  in  both 
X  and  Y  directions,  variable  length  and  size  of  axis  and  plot 
labels,  automatic  or  manual  axis  value  numbering,  horizontal  or 
vertical  orientation,  variable  finished  plot  size  and  the  ability 
to  plot  multiple  curves  and/or  multiple  files  of  data  on  the  same 
plot.  Additionally,  the  module  has  the  ability  to  plot  on  pre¬ 
ruled  graph  paper  by  using  the  reverse  communication  link  to  the 
11T03  for  graph  orientation  data.  A  complete  sample  run  of  the 
Data  Processing  modules  is  contained  in  Section  IV  of  this  report. 

IIl-C.  Program  Development,  Verification  and  Support  Modules 

The  four  modules  in  this  group  provide  general  program 
development  support.  The  first  module  is  the  TI59  optimizer. 

This  module,  as  its  name  implies,  allows  programs  written  for 
the  TI59  to  be  optimized  according  to  known  operational  para¬ 
meters  of  the  calculator.  Designed  originally  as  a  post-processor 
for  the  TI59  Recompiler,  this  module  uses  as  input  a  disk  file 
generated  as  output  of  the  recompiler  or  prepared  with  the  RT-11 
editor  EDIT  or  TECO,  or  human-entered  keyboard  data.  The  opti¬ 
mizer  exploits  known  idiosyncracies  of  the  operation  of  the  TI59 


to  perform  an  "execution- time  program  size"  tradeoff.  This 
time- space  optimization  technique  evaluates  the  tradeoffs  be¬ 
tween  fast  execution  of  in-line  code  versus  the  slower  execution 
of  more  compact  subroutinized  code.  Based  on  operator  input 
parameters,  subroutines  are  forn«d  out  of  repeated  keystroke 
sequences,  literal  digit  strings  are  formed  into  register  based 
constants  where  the  literal  digits  are  replaced  by  a  register 
recall,  and  subroutines  are  organized  in  the  calculator's  memory 
space  according  to  length,  calling  frequency  and  overhead.  The 
resultant  code  executes  as  fast  as  possible  for  a  given  program 
size.  This  module  is  complete. 

To  provide  for  verification  of  translated  programs  and 
assistance  in  the  design  of  Constant  Read  Only  Memories  (CROM) 
for  the  TI59,  a  full  function  simulator  has  been  developed. 

This  module  is,  in  actuality,  a  software  development  aid  in  ad¬ 
dition  to  a  verification  scheme  and  consists  of  a  keyboard  moni¬ 
tor,  multi-mode  debugger,  main  menrary  simulator  and  CROM  memory 
simulator.  The  simulator  provides  all  the  computational  abilities 
of  the  TI59  plus  single  step  execution,  trace  mode,  text  editing 
including  "block  move"  instructions,  a  CROM  paging  scheme  and 
formatted  disk  input/output  including  the  processing  of  manually 
entered  CROM  images.  This  module  is  operational,  however,  several 
enhancements  may  be  desired  to  improve  its  operator  interface  and 
usefulness. 

To  provide  an  automated  basis  for  SIMPLE  program  verifi¬ 
cation,  and  provide  production  translation  capabilities  should 
the  need  arise,  we  have  developed  the  capability  to  translate 
from  SIMPLE  to  RT-11  single  user  BASIC.  Although  not  intended 
as  a  true  MARS  language  processor,  this  module  is  easily  enhanc- 
able  to  that  state.  Input  to  the  program  is  a  MARS  language  pro¬ 
cessor  format  disk  file  containing  SIMPLE  code.  A  file  of  this 
type  is  produced  by  any  of  the  language  decompilers.  The  output 
of  the  processor  is  a  disk  file  that  conforms  to  that  which  is 
expected  by  the  RT-11  BASIC  interpreter.  This  module  is  complete 
and  the  implications  of  its  function  are  obvious.  As  new  genera¬ 
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tions  of  HHPPC's  are  developed  it  is  most  likely  that  many  will 
use  BASIC  rather  than  a  keystroke  based  language.  We  now  have 
the  capability  to  perform  the  necessary  language  translations 
when  the  need  arises. 

The  last  module  in  this  group  is  a  magnetic  card  writing 
aid.  The  hand-printing  of  small  quantities  of  custom  cards  has 
always  been  a  problem.  The  space  allotted  for  labeling  is  small 
so  hand-lettered  labels  become  virtually  unreadable.  This  pro¬ 
gram  uses  the  HP7225A  plotter  and  digitizing  sight,  in  conjunction 
with  the  CRT  to  provide  a  quick  and  convenient  method  of  labeling 
small  quantities  of  magnetic  cards.  This  module  is  complete. 
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This  exsntPle  deoionstrstes  thp  use  of  the  prosrrjms  in  the  cur vef i ttinS 
end  Plotting  PsckasSe#  A  source  sSrsph  from  the  DNA  f'ih- 1  effects  merujel  is 
digitized!*  the  units  sre  trensformed*  s  function^],  form  is  fitted  to  the 
dat3»  end  the  resultant  fit  is  plotted  aSsinst  the  original  for  verifirat; 

Each  step  in  this  examp'le  illustrates  the  capabilities  of  a  different 
program  in  the  peckaSe>  and  also  defiries  a  particular  task  in  the  data 
fitting  process.  The  text  of  this  examr-le  combines  actual  output  from  the 
Programs  with  an  explanation  of  the  results  and  of  the  task  itself.  In 
the  Pioaram  output*  the  user's  responses  to  the  program's  reouests  have 
been  underlined  so  that  they  mas  be  distinguished  from  the  program  materia 
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STEP  1 


Ihe  grsFhicBl  data  is  digitized  and  a  data  file  Froduced  bv 
DAf14  (the  Data  Acouisiton  Module.) 


The  data  to  be  fitted  is  taken  from  the  DNA/EM-I.  effects  manual 
(see  Pis!.  4  >»  Sivind  the  duration  for  the  positive  phase  in  free-air 
overpressure  as  a  function  of  ciistance  from  the/  esr- 1 ‘is i on .  To  convert 
the  drephical  data  into  numeric  data  to  be  ijsed  bv  the  tither  ?  rridrams* 

DAM4  is  used  in  conjunction  with  a  diditizind  tablet. 

The  object  of  riAM4  is  to  produce  a  disk  file  containind  diditized  dat 
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with  DAM4  show  how  this  i 'i  ai  romr-l  isheri! 
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22 


DISTANCE  inwtani 


100  200  500  1000  2000  5000  10000 


DISTANCE  Ifwtl 


Duration  of  Positive  Overpressure  and  Dynamic  Pressure  Phases  for  a  1  kt 
Free  Air  Burst  in  a  Standard  Sea  Level  Atmosphere 


Figure  4. 
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WE  NEED  MINIMUM  AND  MAXIMUM  VALUES  TN  X  AND  Y 

FOR  THE  GRAPH  AS  A  WHOLE >  PLEASE  ENTER  THESE  AS  REAL 

NUMBERS  <WITH  EXPLICIT  DECIMAL  POINT) 

ENTER  XMIN:  100, 

ENTER  XMAX:  10000, 

ENTER  YMIN:  ^ 

ENTER  YMAX;  ^ 

LABEL  FOR  X  AXIS? 

64  CHAR  limit:  DISTANCE  (FT) 

LABEL  FOR  Y  AXIS? 

64  CHAR  limit:  POS,  PHASE  DURATION  (SEC) 

WILL  A  NEXT  DATA  FILE  BE  ASSOCIATED  WITH  THIS  FILE? 

IF  YES»  ENTER  ITS  NAME.  ELSE  ENTER  CARRIAGE  RETURN  ONLY 

FILENAME  ISt  (carriase  return) 

WAS  THERE  A  LAST  DATA  FILE  ASSOCIATED  WITH  THIS  FILE'I’ 

IF  YES.  ENTER  ITS  NAME.  ELSE  ENTER  CARRIAGE  RETURN  ONLY 

FILENAME  IS:  (carriase  return) 

In  order  to  translate  the  internal  coordinates  of  the  disiitizinss 
tablet  into  those  of  the  draeh.  an  algorithm  is  used  which  orients  and 
scales  the  tablet's  output,  based  on  the  designated  coordinater:  of  at 
least  three  points  on  the  tablet,  (Three  points  will  determine  the  avces 
uniauely;  more  than  three  can  be  used  in  an  averasins?  scheme  to  allow  for 
human  imprecision  and  for  noril  inearities  and  skewed  a;:es  in  the  sirar  tis 
themselves , ) 


HOW  MANY  INDIVIDUAL  CURVES  WILL  BE  ON  THIS  GRAPH:  _2 

TO  CALIBRATE  THE  BIT  PAD  WE  WILL  NEED  3  OR  MORE 
WIDELY  SEPARATED  POINTS. -HOW  MANY  SUCH  POINTS  SHALL 
WE  use:  _4_ 

AFTER  EACH  POINT  IS  ENTERED  YOU  MUST  TELL  ME  THE 
DESIRED  X  AND  Y  VALUES  FOR  THf.  POINT, 


INPUT  VALUES  FOR  X  AND  Y:  110, . ,02 


DIGITIZE  A  POINT 


DIGITIZE  A  POINT 
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INPUT  VALUES  FOR  X  AND  Y:  110,».46 


INPUT  VALUES  FOR  X  AND  Yt  920Q ,02 


INPUT  VALUES  FOR  X  AND  YJ  9200*  f ♦ 46 

Now  the  data  values  are  disfitised  directlyi’ 
^?ach  point* 


DIGITIZE  A  POINT 

DIGITIZE  A  POINT 

assisfninsi  a  weight  to 


EACH  CURVE  MIGHT  POSSIBLY  REQUIRE  SEVERAL  PARAMETERS 
FOR  SUBSEQUENT  PROCESSING  BY  OTHER  PROGRAMS.  THE 
PRESENT  PROGRAM  WILL  ACCEPT  WHATEVER  YOU  INPUT... 

THE  SUGGESTED  FORM  IS  ONE  OR  MORE  NUMERICAL  FIELDS 
EACH  ENDING  WITH  A  COMMA  WITH  ANY  DESIRED  CHARACTERS 
FOLLOWING  THE  RIGHTMOST  (OR  ONLY)  COMMA.,... 

PLEASE  ENTER  DATA  CMNTR'OL  FIELD  FOR  THIS  CURVE  WITH 
6A  CHAR  LIMIT!  l.r  ( OVERF-RESSURE ) 

YOU  MAY  NOW  INPUT  POINTS  FROM  THE  DIGITIZER. 

AFTER  EACH  POINT  YOU  WILL  BE  ASKED  FOR  A  WEIGHTING 
FACTOR,  IF  YOU  ARE  MERELY  TESTING  OR  WISH  TO  DISCARD 
THE  POINT,  YOU  SHOULD  GIVE  A  WEIGHT  OF  ZERO.  MOST 
USUALLY  THE  WEIGHT  WILL  BE  ONE.  A  NEGATIVE  WEIGHT 
WILL  BE  ACCEPTED  AS  A  SIGNAL  TO  EITHER  RECALIBRATE: 
(WEI6HT  =  -1)  OR  TO  END  THIS  F'ARTICULAR  CURVE 
(WEIGHT=~2)  (WITHIN  THE  SAME  FAMILY). 


A  CARRIAGE  RETURN  ALONE  GIVES  A  WEIGHT  OF  ONE. 


SELECT  POINT 

ON 

DIGITIZER 

X  IS 

1.03.89 

Y 

WEIGHT? 

TO  DISK 

SELECT  POINT 

ON 

DIGITIZER 

X  IS 

117.68 

Y 

WEIGHT?  :_1_ 

TO  DISK 
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Y  IS  0.34950 


SELECT  POINT  ON  DIGITIZER 

X  IS  8987.1 

WEIGHT?  n 

TO  DISK 

SELECT  POINT  ON  DIGITIZER 

X  IS  454.17  Y  IS  “0.40455  (Sarbaae) 

WEIGHT?  Izl. 

LAST  POINT  DISCARDED »  END  OF  CURVE  SIGNAL  *** 


The  Frocedure  is  repeated  for  the  second  curve  (dwnaniic  pressure). 
After  the  last  point  has  been  diditised.  DAM4  writes  the  completed  file 
of  data  points  in  a  standard  format  for  use  by  the  other  programs. 


SELECT  POINT  ON  DIGITIZER 

X  IS  9635.8  Y  IS 

WEIGHT?  U 

TO  DISK 

SELECT  POINT  ON  DIGITIZER 

X  IS  50  '4.4  Y  IS 

WEIGHT?  :-2 

LAST  POINT  DISCARDED,  END  OF  CURVE  SIGNAL  *** 

DISK  FILE  WRITTEN  AND  CLOSED,  BYE  FOR  NOW 
STOP  --  IN  PROGRAM  DAM4 


The  resultant  disk  file  I 


FILEIDURPP, IDF  DATE : 30- NOV-79  14104 

NAME ’.POSITIVE  PHASE  DURATION*  OVERPRESSURE  S  DYN  PRESSURE 
TYPE  1100 

X(  MIN,  MAX);  1,OOOOOOOE+02,  1 .000()OOOE+04 

Y(MIN,MAX):  O.OOOOOOOE -01 ,  5 . OOOOOOOE-  01 

X  LABELIDISTANCE  (FT) 

Y  LABEL ;P0S.  PHASE  DURATION  (SEC) 

NEXT  file; 

LAST  file; 


0.41904 


--0.35513  (sarbase) 
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SOF  2 

SOC  1 . »  ( OVERPRESSURE ) 
t ,03891 A6E+02. 

1 . 1767580E102» 

1 ,3443927E+02» 

1 ,5292946E+02» 

1 .739536 lE+02» 

1 .9534042E+02» 
2.2219490E+02» 
2,6044104Ei02» 
2.9878143E+02. 
3.3549152E  +  02!. 
3.9502335E+02» 
4.5333554E+02» 

5. 202741 lE+02, 

6. 1280017E+02J 
7.2178143E+02, 
8.4649774E+02. 
9.7575928E+02» 

1 , 164221 9E+03» 

1 .4192300E+03, 

1 ,6643999E+03, 

1 ,9434152E+03. 
2.3185310E+03> 
2.8140056E+03> 
3.4594993E+03» 
4.4018110E+03» 

5. 576481 4E+03» 
6.6231509E+03 j 
7,6986499E+03, 
8.9871709E+03» 

EOC 

SOC  2..  (DYNAMIC  PRESSURE) 
1 .0430090E+02. 

1 ,2546711E+02, 

1 ,4963704E+02» 

1 .7922218E102» 
2.0650964E+02, 
2,5265778E+02» 
2.9992636E+02, 
3,6216669E102» 
4.3739233E+02» 
5.0608652E+02» 
5.8067072E+02» 
6.7789386E+02» 
8.2979028E+02» 
9.7718799E+02» 

1 .2814858E+03» 
1.7096489E+03> 
2.0660352E+03» 
2.5841709E103* 
3,2320801E+03f 
4. 1302988E+03f 
5.2552368E+03» 


1 .6980615E--01 
1.6761231E--01 
1 .6539657E  -01 
1 .6225327E--01 
1 .5629689E~01 
1 .5224892E--01 
1 .4629205E -01 
1 .3650638E-01 
1 .2677549E-01 
1 . 1897544E- 01 
I ,2137309E-01 
1 .3040265E-01 
1 ,4130{:26E-01 
1 .6152833E-01 
1  ,8174B41E--01 
2.0197943E-01 
2, 1756418E~01 
2. 386094 lE-OI 
2,58R2182E"01. 
2.7717680E-01 
2,9179066E-01 
3.0728772E-01 
3.2274100F-01 
3.3534732F-01 
3, 4 78659 7E -01 
3.575ei50E-()l 
3.6370933E-01 
3.6707789E-0I 
3.6949745E-"0  I 


3.7898505E-01 
3.7757584E-01 
3,761S852E-01 
3.7197620E--01 
3.6692446E-01 
3.514()106E-01 
3,3032629E-03 
3.0264145E-01 
2. 833988  7E -0.1 
2.7083200E-01 
2,686()529E-01 
2. 738389 lE-01 
2,8458017E-01 
2.9541999E~()1 
3. 1442809E--01 
3.3245435E-0] 
3,4604251E-01 
3.6235708E-01 
3,7585759E-01 
3.90241  32E-()1 
4.0182194E-01 
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iHHiiliaiiMi 


6,5159072E+03f  4 . 0971&22E“01 f  1. 

7.7386191E+03»  4 . 1396800E-01 »  1. 

9.6358232E+03»  4 . 1903728E~01 f  1. 

EOC 

EOF 
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STEF  2 


The  diiSitised  date  is  reviewed  and  modified  usiniS  proarem  DERM 
(the  Data  Edit  and  Review  Module). 

DERM  reeds  the  file  and  allows  inspection  of  the  data  and  file 
parameters  J 


DERM  version  2  is  ready.  Type  for  information. 

>READ  DURPP.IDF 

Z  1.00000 

Z  2.00000 

53  points  in  data  buffer. 

::  HEADER 

FILE  J  DURPP.IDF 

DATE  :  30-N00-'79 

TIME  :  16504 

NAME  :  POSITIVE  PHASE  DURATION.  OVERPRESSURE  8  DYN  PRESSURE 

TYPE  ;  CARTESIAN 

XLADEL  :  DISTANCE  (FT) 

XAXIS  :  LOGRITHMIC 

XMIN  5  l.OOE+02 

XMAX  :  l,00E+04 

YLABEL  ;  POS.  PHASE  DURATION  (SEC) 

YAXIS  ;  LINEAR 

YMIN  ;  0,00E-01 

YMAX  :  5.00E-01 

NEXT  : 

LAST  ; 


■  LIST 


X 

Y 

Z 

1 

103.892 

0.169806 

1.00000 

rt 

117.676 

0.167612 

1.00000 

3 

134.439 

0.165397 

.1  .00000 

4 

152.929 

0.162254 

1,00000 

5 

173.954 

0,156297 

1.00000 

6 

195.340 

0.152249 

1.00000 

7 

222.195 

0.146292 

1.00000 

8 

260.441 

0.136506 

1.00000 

9 

298.781 

0.126775 

1.00000 

10 

335.492 

0,118975 

1,00000 

11 

395.023 

0.121373 

1.00000 

12 

453.336 

0.130403 

1.00000 
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I 


13 

520.274 

0.141308 

1.00000 

14 

612.800 

0.161528 

1.00000 

15 

721.781 

0.181748 

1 . 00000 

16 

846.498 

0.201979 

1.00000 

17 

975.759 

0.217564 

1.00000 

18 

1164.22 

0.238689 

1.00000 

19 

1419.23 

0.258822 

1.00000 

20 

1664.40 

0.277177 

1.00000 

21 

1943.42 

0.291791 

1.00000 

22 

2318.53 

0.307288 

1.00000 

23 

2814.01 

0.322741 

1.00000 

24 

3459.50 

0.335347 

1.00000 

25 

4401.81 

0.347866 

1.00000 

26 

5576.48 

0,357582 

1 . 00000 

27 

6623.15 

0.363709 

1.00000 

28 

7698.65 

0.367078 

1.00000 

29 

8987.17 

0.369497 

1.00000 

X 

Y 

I 

30 

104.301 

0.378985 

2.00000 

31 

125.467 

0.377576 

2.00000 

32 

149.637 

0.376189 

2.00000 

33 

179.222 

0.371976 

2.00000 

34 

206.510 

0.366924 

2.00000 

35 

252 . 658 

0.351401 

2 . 00000 

36 

299.926 

0.330326 

2.00000 

37 

362.167 

0.302641 

2.00000 

38 

437.392 

0.283399 

2.00000 

39 

506.087 

0.270832 

2 . 00000 

40 

580.671 

0.268605 

2.00000 

41 

677.894 

0.273839 

2.00000 

42 

829.790 

0.284580 

2.00000 

43 

977.188 

0.295420 

2.00000 

44 

1281.49 

0.3:1,4428 

2.00000 

45 

1709.65 

0.332454 

2.00000 

46 

2066.04 

0.346043 

2 . 00000 

47 

2584.17 

0.362357 

2.00000 

48 

3232.08 

0 . 375858 

2.00000 

49 

4130.30 

0.390241 

2.00000 

50 

5255.24 

0.401822 

2 . 00000 

51 

6515.91 

0.409716 

2.00000 

52 

7738.62 

0.413968 

2.00000 

53 

9635.82 

0.419037 

2.00000 

We  will  discard  the  second  curve  (dHnaniic  pressure)  for  now?  and 
convert  the  distance  <X)  from  feet  to  meteors  J 


>DELETE  30- 


>MULT  .304£ 


30 


.>LIST 


X 

Y 

Z 

1 

31.6662 

0.169806 

1 .00000 

o 

35 . 8676 

0.167612 

1.00000 

3 

40.9771 

0.165397 

1.00000 

4 

46.6129 

0.162254 

1.00000 

5 

53.0211 

0.156297 

1 . 00000 

6 

59.5398 

0.152249 

1.00000 

67,7250 

0.146292 

1 .00000 

8 

79.3824 

0.136506 

1 . 00000 

9 

91.0686 

0.126775 

1  .00000 

10 

102.258 

0. 11EJ975 

1.00000 

11 

120.403 

0.121373 

1 .00000 

12 

138.177 

0.130403 

1 .00000 

13 

158.580 

0.141308 

1 .00000 

14 

186.781 

0.1 61528 

1 .00000 

15 

219.999 

0.181748 

1 . 00000 

16 

258.013 

0.201979 

1 . 00000 

17 

297.411 

0.217564 

1.00000 

18 

354.855 

0 . 238689 

1.00000 

19 

432.581 

0 . 258822 

t . 00000 

20 

507 . 309 

0.277177 

1 .00000 

21 

592.353 

0.291791 

1 . 00000 

22' 

706 . 688 

0 . 307288 

1 .00000 

2  3 

857.709 

0.322741 

1 .00000 

'^4 

A«  » 

1054.46 

0.335347 

1 .00000 

•*)  CT 

A-  %.i 

1341.67 

0.347866 

1 . 00000 

26 

1699.71 

0 . 357582 

1 .00000 

27 

2018.74 

0 . 363709 

1.00000 

28 

2346.55 

0 . 367078 

1 . 00000 

29 

2739.29 

0.369497 

1 .00000 

The  header  is  chansed  accordinsfindly >  and  the  niOdi.f:ied  data  is 
t  !  i  e  new  f  i  1  e  ri  a  m  e  li  1)  R  P  F ' .  2  lil "  I 


:  xlabe:l  distance:  (M) 

>XhAX  3000, 

>FILE.  &URPP.21IE' 

>NAHE.  positive:  PHAFSF:  duration,  OtiERFRESSURE 
::  HEADER 

F-ILE  :  DURPP.2EIF 

DATE  ;  30-N0V--7? 

TIME  ;  16104 

NAME  :  POSITIVE  PHASE  DURA] I ON »  OVERPRESSURE 


31 


TYPE 

XLABEL 

XAXIS 

XMIN 

XMAX 

YLABEL 

YAXIS 

YMIN 

YMAX 

NEXT 

LAST 


CARTESIAN 

DISTANCE  (M) 

LOGRITHMIC 

3.00E+01 

3.00E+03 

POS.  PHASE  DURATION  (SEC) 

LINEAR 

O.OOE-01 

5.00E-01 


Finally^  the  modified  datsf  which  is  ready  to  he  fit»  3S  written  into 


a  disk  file  of  the  standard  format J 


::  WRITE 


7  1,00000 

29  Points  written. 


The  resultinsi  file  is* 


FILE;DURPP,2DF  DATE.*30-N0V-79  16. *04 

name:positive  phase  huration,  overpressure 
type: 100 

X(MINfHAX);  3.000000E+01>  3,OOOOOOE+03 

Y(MIN>MAX):  O.OOOOOOE-Ol »  5.000000E-01 

X  label: DISTANCE  (M) 

Y  LABEL  .’POS,  PHASE  DURATION  (SEC) 

NEXT  file: 

LAST  file: 


1,00000 

3.166618E+01, 
3,586758E+01 » 
4,097709E+01 » 
4.661290E+01 > 
5.302106E+01. 
5.953976E+01» 
6.772501E+01» 
7,938243E+01» 
9,106858E+01» 


1 ,698062E-01 
1 .676123E-01 
1 .653966E-  01 
1 ,622S38E-01 
1 .562969E-01 
1  .522489E-01 
1 ,462920E-01 
1 .365064E~()1 
1 .267755E-01 


1 .022578E+02. 

1.189754E 

-01 , 

1 . 

1 .204031E+02» 

1.213731E 

“01 . 

1 . 

1 .381767E+02. 

1 .304027E 

-01 » 

1 . 

1 .585795E+02» 

1 .413()83E 

“01  f 

1  . 

1 .86781 5E+02» 

1 .615283E 

•01  . 

1 . 

2.199990E+02. 

1 .817484E 

-01  » 

1 . 

2.580125E+02!. 

2.019794E- 

“01  f 

1 . 

2.974114E+02t 

2. 175642E 

■01 . 

1 . 

3,E.48549E  +  02. 

2.386894E 

-01  , 

1 . 

4,325813E+02» 

2.588218E- 

•0.1  . 

1 . 

5.073091E+02. 

2.771768E' 

-01 , 

1 . 

5.923530E+02» 

2. 9 179078- 

-01 . 

1 . 

7.066882E+02f 

3 . 072H77f: 

-01 . 

■i 

X  « 

8.577089E+02» 

3.2274101;- 

•01  . 

1 . 

1 .054455E+03» 

3.353473E 

-01 » 

1  - 

1 .341672Ef03. 

3.478660E- 

-01  f 

1 . 

1 .69971 2E+03f 

5.575S15E- 

01. 

1 . 

2.018736E+03» 

3.637094E 

•<■1  ? 

1 . 

2.346548E+03. 

3.670779E 

01  . 

1 . 

2.739290E+03r 

3.694974E- 

-01  . 

1 . 

EGC 

EOF" 
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The  data  is  exdihined  and  a  functional  form  chosen  for  the  fit. 

Then  pro;3ram  FITTER  is  used  to  obtain  the  coefficients  for  which  that 
function  best  approx'imates  the  .diven  data. 

Fig,  5  shows  a  plot  of  the  data  from  data  file  DURPF' .  I'liF .  (This 
figure  was  produced  by  the  program  HP7225.) 

A  functional  form  which  suggests  itself  is 
a  d 

Y  =  - - -  +  - - - 

(X/b)  +  1  (X/g)~"  +  1 

where  Y  =  positive  phase  duration  in  seconds 
X  -  distance  from  burst  in  meters. 

The  term  a  /  C(X/b)^  +  13  governs  the  limit  of  small  X  (X  <<  brg) 
and  defines  the  left  side  of  the  graph?  likewise  the  term 

d  /  r(X/g)  +  13  governs  the  limit  of  large  X  (X  »  t'»g)  and  defines  the 
right  side  of  the  graph,  (See  fig,  5.) 

A  1 inearisation  of  parameters  technioue  is  used  to  determine  the 
coefficients  at  b>  cr  d»  g>  and  h.  This  method  reeuires  initial  estimates 
for  the  values  of  the  coeff icients»  which  we  may  get  by  cvsamining  the 
asymptotic  behavior  of  the  function. 

For  X  <<  b»g»  Y  tends  to  a  ™  ,17  seconds  (see  fig,  4  or  5) 

For  X  >>  b»g>  Y  tends  to  d  =  ,375  seconds 

When  X  =  bf  Y  =  a/2  (neglecting  the  second  term  in  the  function) 

which  we  estimate  to  occur  at  X  ~  b  =  400  ft  =  122  meters. 

When  X  =  g»  Y  =•  d/2  (neglecting  the  first  term  in  the  function) 

which  we  estimate  to  occur  at  X  =  g  =  800  ft  =  244  meters, 

Finallyf  as  initial  guesses  for  coefficients  c  and  h  we  cJioose  1, 
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DISTANCE  (M> 


Now  F'roSirsni  FITTER  is  used  to  optimize  these  coef  f  icients  *  Firsts 
the  functional  form  to  be  used»  which  is  rather  seec  i al.  i zed ?  is  set  ue 
in  a  FORTRAN  subroutine  to  be  called  by  the  main  erosram#  This  subroutine 
will  be  used  for  each  of  the  two  terms  in  the  furiction. 


FUNCTION  USER3(l!.X)  !  UBFP?. .  RPD 

C 

C  USER-liEFINED  F^OUTINE  FOR  POS.  PHASE  DURATION  FIT 

C 

DOUBLE  PRECISION  PI 

REAL  COEFF  (  4  » 1 3  >  1 0  )  i-  DATA  ( 2  -  4  j  1 01  ) ,  SCRTCH  (  3 ,2,1  S  i 
INTEGER  FUNCT(4, 14) , 1NF0( 10) 

BYTE  CLEAR ( 7 ) , CHAR , CHR ( 4  > , JUNK 
REAL  EXT(6) 

C 

COMMON/ALL/COEFF  ,  DATA ,  SCRTCH ,  FUNCT  ,  I NFO  ,  CLEAR  ,  CHR  ,  JUNK  ,  P I  ,  ZPC)  I  NT 
COMMON/USER/EXT 
C 

DATA  EXT(3)/'PPD'/ 

C 

C  F(X)  -  A  /  ((X/B)"'C  T  1) 

C 

USER3--C0EFF  (1,1,1  > ./  <  ( X/COEFF  (.1,1,2))  **C0EFF  (1,1, 3  )  f  1 .  ) 

C 

RETURN 

END 


The  procedure 

used  to 

opt i miz 

e  the 

coetf  icier  its  is 

an  itcM'at 

i ve  one 

which 

reouires  that 

some  of 

the?  coefficients  he 

held  Pi 

xed  while 

the  other 

are  allowed 

to  vary 

.  If  St 

f  i  rst 

t h e  c  o  e  f  f  :i  c  i  e  r vt  <■.  a  ,  1.' 

,  d ,  and  S 

are  fixe 

and  c 

and  h 

are  optimized. 

we  Set 

the  followiriS 

resu 1 ts 

fi  om  tlie 

fitter: 

iteration 

a 

b 

c 

d 

S 

h 

RMS  rels 

Live  e  r  r-  o 

0 

,17 

1 22 , 

1.000 

.375 

244. 

1 ,000 

1 

.17 

122. 

3,319 

.375 

244, 

2, 199 

,1313 

o 

,  17 

1 22 , 

3.022 

.375 

244 . 

2,155 

,1307 

3 

,17 

122. 

2,830 

.375 

244, 

2,125 

,  1  309 

4 

.17 

122, 

2.935 

.375 

244, 

2,137 

,  1315 

5 

.17 

122. 

2.879 

,375 

244, 

2,130 

,  1315 

Now,  fixind  the  values  of  c  and  h  as  they  have  stabilized,  and 


allowing  3»  b»  d»  and  s  to  vary 


iteretion 

a 

b 

c 

d 

£i 

li 

RMS  relative  erro 

6 

.  1641 

104.4 

2.879 

.  3565 

253 . 1 

2.130 

.0421 

7 

.1675 

101.6 

2 . 879 

.3564 

253 . 2 

2.  130 

.0409 

8 

.1679 

101.2 

2.079 

.  3564 

252.8 

2.130 

.0407 

9 

.1679 

101.1 

2.879 

.  3563 

252 . 7 

2.130 

.0406 

Finally? 

c  and 

h  are 

as^ain  allowed 

to  vary? 

to  achieve  the 

optinium  vsluesJ 


iteration 

a 

b 

c 

d 

S 

h  RMS 

relative  erro 

10 

.1641 

104.4 

2.552 

.3565 

253 . 

1  1.978 

.0415 

11 

.  1 64 1 

104.4 

2.594 

.  3565 

V-  5  3  * 

1  1.991 

.0406 

12 

.  1641 

104.4 

2 . 590 

.  3565 

253 . 

1  1.990 

.0406 

with  stable 

values . 

Prosra 

ir.  FITTER 

Si  ves 

the 

PD  int--by -point 

error 

analysis  5 

I 

X 

Y 

F 

F-Y 

(F-Y)/Y 

1 

3.16662E+01 

1 .69806E--01 

1.65635E-01 

-4. 17146E-03 

-2.45660E- 02 

o 

3.58676E+01 

1 .67612E-01 

1  ..64364F-01 

-"3.24813E-03 

-1 .93788F-02 

3 

4.09771E+01 

1 .65397E-01 

1 .62452E-01 

-2.94496E-03 

-  1 .78055E-02 

4 

4.66129E+01 

1 .62254E-01 

1 .5992  IE -01 

-2.33309E-03 

-1 .43793E-02 

5 

5.30211E+01 

1 .56297E-01 

1 .56612E-01 

3. 15145E-04 

2.01632E-03 

6 

5.95398E+01 

1 .52249E-01 

1 .52937E-01 

6.87972E-04 

4.51873E-03 

7 

6.77250E+01 

1 .46292E-01 

1 .48172E-01 

1 .88008E-03 

1 .28516E-02 

8 

7.93824E+01 

1 .36506E-01 

1 .41775E-01 

5.26822E-03 

3.85932E-02 

9 

9. 10686E+01 

1 .26776E-01 

1 .36586E-01 

9.81008E  -03 

7.73815E  -02 

10 

1 .02258E+02 

1 . 18975E-01 

1 .33253E-01 

1 .42777E-02 

1 .20005E  -01 

11 

1 .20403E+02 

1 .21373E-01 

1 .31606E-01 

1 .02327E-02 

8.43075E-02 

12 

1 .38177E+02 

1 .30403E-01 

1 .34130E-01 

3.72693E-03 

2.85802E-02 

13 

1 .58579E+02 

1.41308E-01 

1 .40914E-01 

-3.94389E-04 

-2.79098E-03 

14 

1 .86781E+02 

1 .61528E-01 

1 .54592E-01 

-6.93622E-03 

-4.29412E--02 

15 

2. 19999E+02 

1 .81748E-01 

1 .73536E-01 

-8.21231E-03 

-4.51851E- 02 

16 

2.58012E+02 

2.01979E-01 

1 .95495E-01 

•-6.484e7E-03 

-3.21066E  02 

17 

2.97411E+02 

2. 17564E-01 

2. 16485E-01 

-1 .07932E  -03 

-4.96093E-03 

18 

3.54855E+02 

2.38689E-01 

2.42435E-01 

3.74570E-03 

1 .56928E-02 

19 

4.32581E+02 

2.58822E-01 

2.6912RE-01 

1 .03061E- 02 

3.98192E-02 

20 

5.07309E+02 

2.77177E-01 

2.87652E--01 

1 ,04756E-02 

3.77939E-02 

21 

5.92353E+02 

2.91791E-01 

3.02795E~01 

1 . 10040E-02 

3.77120E--02 

22 

7.06688E+02 

3.07288E-01 

3. 16665E-01 

9.37760E--03 

3.05173E-02 
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23 

8.57709E+02 

3.22741E-01 

3.28223E-01 

5.48199E-03 

1 .69857E-02 

24 

1 .05446E+03 

3.35347E-01 

3.371 14E-0i 

1.76695E-03 

5.26901E  -03 

25 

1.34167E+03 

3.47866E-01 

3.44144E-01 

-3.72222E-03 

-1 .07002E-02 

26 

1 .69971E+03 

3.57581E-01 

3.48603E-01 

-■8.97e81E-03 

-2.51098E-02 

27 

2.01874E+03 

3.63709E-01 

3.50801E-01 

•  1 .2908  IE -02 

-3.54901E- 02 

28 

2.34655E+03 

3.67078E-01 

3.52212E--01 

-1  .48655E--02 

-4 .04967E-02 

29 

2.73929E+03 

3.69497E-01 

3.53293E-01 

-1 .62042E-02 

•4.38548E  -02 

ERROR  is:  4.06593E-02  RMS  REL  ERROR 


Since  this  eaustion  will  be  used  in  a  calculator  prodram  with  limited 
space  available;  it  is  important  to  minimize  the  number  of  significant 
digits  reeuired  in  the  coefficients.  The  FITTER  allows  us  to  round  the 
values  of  the  coefficients;  and  examine  the  resulting  errors.. 


a 

b 

c 

d 

g 

h 

.  17 

100. 

2.6 

.36 

260. 

— ■  ^ 

ith  the  errors: 

1 

X 

Y 

F 

F-Y 

(F-Y)/Y 

1 

3. 16662Ef01 

1 .69806E-01 

1 .67121E-01 

-2.6B531E-03 

-1 .58140E-02 

o 

3.58676E+01 

1 .67612E-01 

1 .65670E-01 

-1 .94208E  -03 

-1 . 15868E-02 

3 

4.09771E+01 

1 .65397E-01 

t.6350aE-01 

-1 .8882BF-03 

-1 . 14167E -02 

4 

4.66129E+01 

1 .62254E-01 

1 .60669E-01 

-1 .5B4R9E-03 

-9. 76795E-03 

5 

5.30211E+01 

1 .56297E-01 

1 .56977E-01 

6.79B80E-04 

4.34993E-0;3 

6 

5.95398E+01 

1 .52249E-01 

1 .52889E--0] 

6.40094E-04 

4.20426E-03 

7 

6. 77250E+01 

1 .46292E-01 

1 .47596E-0] 

) ,30381F-03 

e.91236E- 03 

8 

7.93824E+01 

1 .36506E-  01 

1.40471E-01 

3 . 964B0l.-:-03 

2.90448E-02 

9 

9. 10686E+01 

1 .26776E-01 

1 .34627E-01 

7.85197E-03 

6. 19360E--02 

10 

1 .02258E+02 

1  . 18975E-01 

1 .30760E-01 

1 . 17846E-02 

9.9051  IE- 02 

11 

1 ,20403E+02 

1 .21373E-01 

1 .28442E-01 

7,06901 E -03 

5.8242OE-02 

12 

1 .38177E+02 

1 .30403E-01 

1 .30519E-01 

1 . 16229E-04 

8.91309E-04 

13 

1 .58579E+02 

1 .41308E-01 

1 .36996E-01 

-4.31,1 84E-03 

-3.05137E-02 

14 

1 .86781E+02 

1 .61528E-01 

1 .50528E-01 

-  1 . 10002E -02 

-6.81009E-02 

15 

2. 19999E+02 

1 ,81748r-01 

1.69595E-01 

-  1 .21530E-02 

-6.6a672E-02 

16 

2.58012E+02 

2.01979E-01 

1 .919461: -01 

-1 .00337F-02 

-  4.96771E-02 

17 

2.97411E+02 

2. 17564E-01 

2. 13492F-01 

-4.0719BE  -03 

-1 .87162E-02 

18 

3.54855E+02 

2.38689E~01 

2.40335E-  01 

1 .64548E  -03 

6.89379E:  -03 

19 

4.32581E+02 

2.58822E-01 

2.68153E-01 

9.331 43F-n3 

3.60535E -02 

20 

5.07309E+02 

2,77177E-01 

2.8756HE-01 

1  .0391 7F  -02 

3.7491  2F.  -02 

21 

5.92353E+02 

2.91791E-01 

3.03497F-01 

1 . 17064E-02 

4.011 90F-02 

02 

7.06688E+02 

3.07288E-01 

3. 18127E-01 

1 .083M9E-02 

3.527281.  02 

23 

8.57709E+02 

3.22741E-01 

3.3033RF-01 

7.596  70E  -03 

2.353818-02 

24 

1 .05446E+03 

3,35347E-01 

3.39738E-01 

4 . 391 04F -03 

1  .,J0940E-02 

25 

1 .34167E  +  03 

3.47866E-01 

3.47169E-01 

-6.97374F-04 

2.004  '2r- 03 

26 

1 .69971E+03 

3.57581E-01 

3.51876F- 01 

-5. 70S06F  0  1 

-1.59f,4AF  02 

27 

2.01874E+03 

3.63709E-01 

3.54195F  01 

9. 5 148 IF  -o  < 

2 .616058-02 
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28  2.34655E+03  3.67078E-01  3.55680E-01  -1.13975E-02  -3.10493E-0? 

29  2«73929E+03  3.69497E-01  3.56817E-01  -l,2680AE-02  -3.43184E~02 

ERROR  is;  3.7547('.E-02  RMS  REL  ERROR 

which  are  General ly  acce^-tahle* 


STEP  A 


A  Plot  of  the  fitting  function  is  made  sdeinst  the  original  graph 
to  visually  verify  the  fit« 

Using  programs  FITTER  and  DERM*  a  data  file  containing  the  values 
of  our  fitting  function  over  the  range  of  the  source  graph  was  crestedJ 


FILEtDATA. IDF  DATE : 05-DEC-79  13:55 

name:pos.  phase  fit*  overpressure 
type: 100 

X<MIN*MAX):  1 .OOOOOOE+02*  l.OOOOOOE+04 

Y(MIN*MAX):  O.OOOOOOE-01*  5.000000E-01 

X  label: DISTANCE  (FT) 

Y  LABEL :P0S.  PHASE  DURATION  (SEC) 

NEXT  file: 

LAST  file: 


SOF 

SOC  0.000000 

1 .OOOOOOE+02* 
1 ,047615E+02* 
1 .097498E+02. 
1.149756E+02* 
1.204502E+02* 
1 .261855E+02* 
1 .321939E+02* 
1.384883E+02* 
1.450825E+02* 


6.579140E+03> 

6.892408E+03* 

7.220593E+03* 

7.564404E+03* 

7.924587E+03* 

8.301918E+03* 

8.697218E+03* 

9.111340E+03* 

9,545179E+03* 

9.999677E+03* 

EOC 

EOF 


1 ,674751E-01 *  1. 
1.A70384E--01*  1. 
1.665415E-01*  t, 
1.659782E-01*  1. 
1 ,653418E-01 *  1. 
1 .646254E-01 *  1. 
1.638220E-01*  1. 
1.629248E-01 >  1. 
1 ,619272E-01 *  1. 


3.541182E-01 

3,546310E-01 

3.550998E-01 

3,555282E-01 

3.559196E-()1 

3.5A2773E--01 

3.5/.6038E-01 

3.569020E“01 

3.571743E-()1 

3.574228E-01 


1, 
1 . 
1 . 
1 . 
1 . 
1 . 
1  . 
1  . 
1 . 
1. 


The  plotter  program  HP7225  scales  and  orients  the  axes  for  the  plot 


40 


m  accordance  with  the  data  file  esrameters  and  additional  user  input 


ENTER  PLOT  FILE  NAME*.  DATA .  1 DF 
OPTIONS  ?  ?_ 

THE  FOLLOWING  OPTIONS  ARE  AVAILABLE 


SYM  -----  -  SYMBOL  l-T.  PLOT  MODE 

ROT  -  --  - .  90  DEGREE  PLOT  ROTATION 

LINE - CONTINUOUS  LINE  MODf: 

AXES - -  CHANGE  AXES  TYPES 

SIZE  -----  -  CHANGE  GRAPH  SIZE 

MIN - CHANGE  GRAPH  MIN  OR  MAX 

MAX  -----  CHANGE  GRAPH  MIN  OR  MAX 

INFO - -  -  -  DISPLAYS  F'ROGRAM  USER  INFORMATION 

SET - -  -  SET  GRAPH  AXES  FRAME  LIMITS 


OPTIONS  ?  ROT 
OPTIONS  ?  S£r 

PLEASE  POSITION  PEN  FOR  GRAPH  LOWER  LEFT  CORNER- 
HIT  •RETURN'  OR  'ENTER'  WHEN  POSITIONED 

PEN  position: 

X  =  4626 

Y  =  7000 

PLEASE  POSITION  PEN  FOR  GRAPH  UPPER  RIGHT  CORNER- 
HIT  'RETURN'  OR  'ENTER'  WHEN  POSITIONED 

PEN  position: 

X  =  7346 

Y  =  1520 


The  curve  itself  is  now  Plotted^  with  a  dotted  line: 


OPTIONS  ?  (carriasJe  rectum) 

READY  FOR  CURVE  ♦  1  ?  _l^  (indicates  a  reouest  for  l  ine  tyee) 

ENTER  LINE  TYPE  0  THROUGH  6, 

FOR  CONTINUOUS  LINE.  JUST  HIT  A  CARRIAGE  RETURN 
LINE  TYPE?  ===>  0 


The  curve  is  tJu-.i  plotted?  and  the  result  is  fis.  6.  A  coniearison 
of  the  orisfinal  digitized  data  with  the  final  fit  maw  also  be  made?  with 
the  result  of  fiH,  7  (the  solid  line  represents  the  fitting  function?  whil 
the  dots  are  individual  digitization  points).  We  have  arrived  at  a  fairlw 
accurate  fit?  with  pleasin^lw  short  coefficients. 


41 


POSITIVE  PHASE  DURATION  itite) 


DISTANCE  (nwwit) 


100  200  500  1000  2000  5000  10000 

DISTANCE 

Duration  of  Positive  Overpressure  and  Dynamic  Pressure  Phases  for  a  1  kt 
Free  Air  Burst  in  a  Standard  Sea  Level  Atmosphere 


Figure  6 
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Studied  &  Analyses 

Department  of  Energy 

Department  of  the  Air  Force 

ATTN:  CTID 

ATTN 

AF/SAG,  H.  Zwemer 

ATTN 

AF/SD,  C.  Cupp 

DEPARTMENT  OF  ENERGY  CONTRACTORS 

ATTN 

AF/SAS 

Lawrence  Livermore  National  Lab 

Deputy  Chief  of  Staff 

ATTN:  L-IO,  H.  Kruger 

Operations  Plans  A  Readiness 

AHN:  L-96 

Department  of  the  Air  Force 

ATTN:  L-49,  W.  Hogan 

AHN 

AFXOOIR 

ATTN;  L-24 

ATTN 

AFXOXFS 

ATTN:  L-21 

ATTN 

AFXOXFM 

ATTN:  L-9,  R.  Barker 

ATTN:  Technical  Info  Dept  Library 

Air  Force  Intelligence  Service 

ATTN;  L-153,  H.  Hofer 

ATTN;  I  NOT 

ATTN:  L-122,  D.  Manner 
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DEPARTMENT  OF  ENERGY  CONTRACTORS  (Continued) 

Los  Alamos  National  Lab 

ATTN:  MS  670.  J.  Hopkins 

ATTN;  MS  218,  R.  Seldon 

ATTN;  MS  670,  J.  McQueen 

ATTN:  R.  Selden 

ATTN:  MS  948.  A.  Charmatz 

ATTN:  Reports  Library 

ATTN:  Sandoval /Chapin/Lyons/Best/Dowler 

ATTN:  MS  634,  T.  Dowler 

Lovelace  Biomedical  i  Environmental  Rsch  Inst,  Inc 
ATTN:  D.  Richmond 

Sandia  National  Labs 

Livermore  Lab 


DEPARTMENT  OF  DEFENSE  CONTRACTORS  (Continued) 

Effects  Technology,  Inc 
ATTN :  M.  Rosen 

Ford  Aerospace  &  Communications  Corp 
ATTN:  K.  Attinger 

General  Electric  Co 

ATTN:  J.  Hannabach 

General  Research  Corp 
ATTN:  8.  Gragg 

General  Research  Corp 

ATTN;  Tactical  Warfare  Operations 


ATTN: 

8324,  J.  Strauve 

Horizons  Technology,  Inc 

ATTN: 

J.  Wirth 

ATTN 

A.  Bemath 

ATTN: 

H.  Norris 

ATTN 

R.  Kruger 

ATTN 

J.  Shaibly 

Sandia  National  Lab 

4  cy  ATTN 

W.  Caid 

ATTN: 

5613,  R.  Stratton 

4  cy  ATTN 

W.  Simmons 

AHN: 

A.  Chabai 

ATTN: 

C.  Mehl 

Hudson  Institute.  Inc 

ATTN; 

3141 

ATTN 

H.  Kahn 

ATTN; 

J .  PI  impton 

IIT  Research  Institute 

DEPARTMENT  OF  DEFENSE  CONTRACTORS 

ATTN 

Documents  Library 

Advanced  Research  i  Applications  Corp 

Institute  for  Defense  Analyses 

ATTN; 

Document  Control 

ATTN 

J.  Grote 

ATTN 

Classified  Library 

Aerospace  Corp 

ATTN: 

R.  Crolius 

IRT  Corp 

ATTN: 

Library 

ATTN 

N.  Rudie 

Analytical 

Assessments  Corp 

JAYCOR 

ATTN; 

A.  Wagner 

ATTN 

R.  Stahl 

Aptek,  Inc 

JAYCOR 

ATTN: 

T.  Meagher 

ATTN 

R.  Sullivan 

AVCO  Research  &  Systems  Group 

JAYCOR 

ATTN: 

P.  Grady 

ATTN 

R.  Schaefer 

BDM  Corp 

Johns  Hopkins  University 

ATTN: 

J.  Conant 

ATTN 

F.  Dreisford 

ATTN: 

J.  Braddock 

ATTN: 

J.  Bode 

Kaman  Avidyne 

ATTN 

N.  Hobbs 

Beers  Associates,  Inc 

ATTN: 

B.  Beers 

Kaman  Sciences  Corp 

ATTN 

D.  Seitz 

Boeing  Co 

ATTN 

F.  Shelton 

ATTN: 

B.  Lempriere 

ATTN 

F.  Rich 

Boeing  Co 

Kaman  Tempt 

AHN: 

0.  Choate 

ATTN 

DAS  1  AC 

Commander 

KTECH  Corp 

66th  MI  Group 

ATTN 

D.  Keller 

2  cy  ATTN: 

T.  Greene 

DEPARTMENT  OF  DEFENSE  CONTRACTORS  (Continued) 

Lockheed  Missiles  &  Space  Co,  Inc 
ATTN:  R.  Walz 
ATTN:  J.  Fortun 

Martin  Marietta  Corp 
ATTN:  F.  Marion 


DEPARTMENT  OF  DEFENSE  CONTRACTORS  (Continued) 

Santa  Fe  Corp 

ATTN:  D.  Paolucci 

Science  &  Engrg  Associates,  Inc 
ATTN:  J.  Cramer 

Science  Applications,  Inc 


Martin  Marietta  Corp 

ATTN 

L.  Nesseler 

ATTN 

J.  Donathan 

ATTN 

W.  Yengst 

ATTN 

Document  Control 

Maxwell  Labs,  Inc 

ATTN 

J.  Warner 

AHN 

J.  Rauch 

ATTN 

M.  Drake 

ATTN 

U.  Clark 

Science  Applications,  Inc 

McDonnell  Douglas  Corp 

ATTN 

J.  Foster 

AnN 

L .  Cohen 

ATTN 

W.  Layson 

AHN 

D.  Newlander 

ATTN 

B.  Bennett 

ATTN 

R.  Craver 

Mission  Research  Corp 

ATTN 

Doc  Control 

ATTN 

J.  Gilbert 

ATTN 

J.  Raymond 

Science  Applications,  Inc 

ATTN 

D.  Kaul 

Mission  Research  Corp 

ATTN 

V.  Van  Lint 

SRI  International 

ATTN 

D.  Curran 

Pacific-Sierra  Research  Corp 

ATTN 

P.  Dolan 

ATTN 

H.  Grade 

SRI  International 

PDA  Engineering 

AHN 

J.  Scholz 

ATTN 

P.  Gilliam 

ATTN 

J.  SI  OSS 

Physical  Sciences,  Inc 

System  Planning  Corp 

ATTN:  P.  Spence 

ATTN 

W.  Robertson 

Physics  International  Co 

Systems,  Science  &  Software,  Inc 

ATTN 

P.  Spangler 

ATTN 

K.  Pyatt 

ATTN 

J.  Shea 

ATTN 

C.  Dismukes 

R  &  D  Associates 

T.N.  Dupuy  Associates,  Inc 

ATTN 

S.  Borjon 

ATTN 

T.  Dupuy 

ATTN 

R.  Port 

ATTN 

A.  Field 

Tetra  Tech 

Inc 

AnN 

W.  Graham 

ATTN 

F.  Bothwell 

ATTN 

C.  Rogers 

ATTN 

P.  Haas 

TRW  Electronics  &  Defense  Sector 

ATTN 

A.  Lynn 

ATTN 

H.  King 

ATTN 

P.  Rausch 

ATTN 

W.  Polich 

ATTN 

E.  Carson 

TRW  Electronics  &  Defense  Sector 

Rockwell  International  Corp 

ATTN 

R.  Burnett 

ATTN 

G.  Messenger 

ATTN 

D.  Scally 

ATTN 

N.  Lipner 

AnN 

D.  Clement 
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